Molecular correlations and solvation in simple fluids 
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We study the molecular correlations in a lattice model of a solution of a low-solubility solute, 
with emphasis on how the thermodynamics is reflected in the correlation functions. The model is 
treated in Bethe-Guggenheim approximation, which is exact on a Bethe lattice (Cayley tree). The 
solution properties are obtained in the limit of infinite dilution of the solute. With hn(r), /112(f), 
and /122(f) the three pair correlation functions as functions of the separation r (subscripts 1 and 2 
referring to solvent and solute, respectively), we find for r > 2 lattice steps that Zi22(f)//ii2(f) = 
/ii2(f)//in(f )■ This illustrates a general theorem that holds in the asymptotic limit of infinite r. 
The three correlation functions share a common exponential decay length (correlation length), but 
when the solubility of the solute is low the amplitude of the decay of /122(f) is much greater than that 
of /112(f), which in turn is much greater than that of /in(f). As a consequence the amplitude of the 
decay of /122(f) is enormously greater than that of /in(f). The effective solute-solute attraction then 
remains discernible at distances at which the solvent molecules are essentially no longer correlated, 
as found in similar circumstances in an earlier model. The second osmotic virial coefficient is large 
and negative, as expected. We find that the solvent-mediated part W(r) of the potential of mean 
force between solutes, evaluated at contact, r = 1, is related in this model to the Gibbs free energy of 
solvation at fixed pressure, AGp, by (Z/2)W (1) + AGp = pvo, where Z is the coordination number 
of the lattice, p the pressure, and vo the volume of the cell associated with each lattice site. A large, 
positive AGp associated with the low solubility is thus reflected in a strong attraction (large negative 
W at contact), which is the major contributor to the second osmotic virial coefficient. In this model, 
the low solubility (large positive AGp) is due partly to an unfavorable enthalpy of solvation and 
partly to an unfavorable solvation entropy, unlike in the hydrophobic effect, where the enthalpy of 
solvation itself favors high solubility, but is overweighed by the unfavorable solvation entropy. 

PACS numbers: 



I. INTRODUCTION 

When a solute is only sparingly soluble in a solvent, 
pairs of solute molecules effectively attract. In the ex- 
treme, the solutes aggregate and precipitate as a separate 
phase in equilibrium with a very dilute solution of that 
solute in the solvent. We propose here to obtain quanti- 
tative relations between the solute's low innate solubility 
and the potential of mean force between a pair of so- 
lute molecules, the latter describing the effective attrac- 
tion between them. Of particular interest is the solvent- 
mediated part of the potential of mean force; i.e., the 
part remaining after the direct solute-solute interaction 
is subtracted. 

We study this in a lattice model on a Bethe lat- 
tice (Cayley tree), for which the Bethe-Guggenheim 
"quasichemical" approximation 0, Q becomes exact. 
The Bethe-Guggenheim and related approximations have 
been used before for lattice models of water and hy- 
drophobic solvation and the approximation is used 
extensively in continuum models of solutions as well [9, 
Hoj . Models of water and aqueous solutions on other 
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hierarchical lattices have also been extensively stud- 
ied [ill, (HI- The equivalent of the present model but 
with emphasis on its properties in the neighborhood of 
its critical-point manifold rather than in the thermody- 
namic states of interest here has also been analyzed pre- 
viously fl3j . Here we are primarily interested in the limit 
of infinite dilution of the solute in the dense liquid sol- 
vent. 



In the following section we review and collect the prop- 
erties of the pure solvent on the Bethe lattice, which we 
will need subsequently. In Section IIHI we find the ther- 
modynamic properties and correlation functions of the 
solution in the limit of an infinitely dilute solute. These 
are the main results in the paper, and are illustrated nu- 
merically with a representative set of the model's param- 
eters. Section IIVI re-expresses the results in two limits, 
both helping to illuminate the general formulas. One is 
the mean-field limit, in which the coordination number 
Z of the lattice becomes infinite while the solvent-solvent 
interaction energy —en goes to 0, with Zen fixed. The 
other is the low-temperature limit, T — > 0, while the 
high-density liquid solvent remains in coexistence with 
its low-density equilibrium vapor. The final Section IVl is 
a brief summary of our results. 
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II. SOLVENT PROPERTIES 

Here we summarize the properties of the one- 
component system, without solute, which will then be 
the properties of the pure solvent we shall need in the 
later parts of the paper. The formulas here are either 
well known or readily derivable and they establish our 
notation. 

Our Bethe lattice is of coordination number Z. As- 
sociated with each site is a cell of volume Vq [l3J ■ Only 
singly occupied and empty cells are allowed. Except for 
the exclusion of multiple occupancy, the only interactions 
are between molecules in neighboring cells, which inter- 
act with energy —en (en > 0). With p the density, T 
the temperature, p the pressure, and k Boltzmann's con- 
stant, we introduce the abbreviations 

y = pv , t = kT/en, TT=pv /en, c = . (1) 

The 7r, y, t equation of state on the Bethe lattice ("qua- 
sichemical" approximation) is then 



t 



where 



with 



_ 2y-l + Q 
a ~ 2(1 -y)c 

= y/l + 4(c-l)(l-y)y. 



(2) 



(3) 



(4) 



We define a dimensionless compressibility \ an d coef- 
ficient of thermal expansion 9 by 
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(5) 



these being, respectively, the fractional rate of increase of 
density with increasing pressure at fixed temperature and 
the fractional rate of decrease of density with increasing 
temperature at fixed pressure. In the high-density states 
of interest here, these will both be very small (<C 1), as 
in real liquids. 

The equation of the liquid-vapor coexistence curve in 
the y,t plane is 
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The equation of the spinodal is 
Z 



Q 



Z-2 



(spinodal). 



(coex. ). 

(6) 

(7) 



The spinodal lies within the coexistence curve in the y,t 
plane and is tangent to it at the critical point. The crit- 
ical point j/c, t c is at 
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(8) 



We measure distances r between cells as the number of 
steps on the underlying lattice. This defines it uniquely 
because there are no closed loops on a Bethe lattice (Cay- 
ley tree). The pair correlation function hu(r) is 

h n (r) = (l-ije-** (r>l), /i u (0) = -l, (9) 

where the correlation length (exponential decay length) 
£is 
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1 



(10) 



When, as in © and (|TU|) with Q as in the thermo- 
dynamic state is specified by the variables y and t, the 
correlation function is independent of coordination num- 
ber Z and so is the same as when Z — 2 (the linear 
chain). That is a reflection of the uniqueness of the path 
between any two sites on the lattice. But it is no longer 
true when some other variable such as 7r in Eq. ([2]) , or the 
activity, replaces either or both of y and t] then h\\{r) 
does depend on Z . Likewise, when the state is specified 
to be at the coexistence curve, hn(r) again depends on 
Z via ([6]). Indeed, there is no coexistence curve unless 
Z > 2. 

It is a peculiarity of the Bethe lattice that the correla- 
tion length £ in (|10p remains finite at the critical point. 
Nevertheless, the sum of hn(r) over all lattice sites di- 
verges there, along with the compressibility, in accord 
with the Ornstein-Zernike theorem. That is because the 
number N r of sites distant r from any central site diverges 
exponentially rapidly as r — > oo: 

N r = Z(Z-l) r - 1 (r>l), JV = 1. (11) 

Thus, as may be verified from ©-© and (|8))- (|TTt . 



t X /n = l/y + J2 N rhu(r), 



(12) 



r=0 



as required by the Ornstein-Zernike theorem. Both sides 
diverge at the critical point given by ©. 



III. SOLUTION WITH INFINITELY DILUTE 
SOLUTE 

We define the solubility E of the solute (species 2) to 
be the dimensionless ratio of its number density p2 in the 
solution to its activity in the limit of infinite dilution, 
with the activity defined to become equal to its number 
density in an ideal gas; thus, 



P2/z 2 - 



(13) 



Except for a factor of kT, this is also the Henry's-law 
constant, and it is also the Ostwald absorption coefficient 
when the vapor in equilibrium with the solution is an 
ideal gas. Low solubility means E 1. 
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By the potential-distribution theorem (particle- 
insertion theorem) , E is then the product of the probabil- 
ity 1—y that a trial solute molecule inserted into the pure 
solvent has found an empty cell, multiplying the condi- 
tional average (e~*/ fcT ), with ^ the energy of interac- 
tion of the inserted solute with its (hypothetically undis- 
turbed) neighboring solvent molecules. On the Bethe 
lattice of coordination number the probability fj that 
an empty cell in the pure solvent is neighbored by exactly 
j occupied cells is 



Z\ 



j\{Z-j)\{l + a) z ' 
with a as in ([3]). Thus, 

E = (l-zy)X>e jei2/fcT 

z 



(14) 



(1-1/) 



1 + ac a 
1 + a 



(15) 



where — e\2 is the interaction energy between solute and 
solvent molecules in neighboring cells, where c is as in ([1} , 
and where 



ei2/e 



ii- 



(16) 



This gives E as a function of temperature t and solvent 
density y. In the dense liquid solvent at given t it is a 
very sensitive function of y but insensitive to pressure 
7r because of the solvent's low compressibility. We shall 
later wish to know the rate at which E changes with 
temperature at fixed pressure, which may be obtained 
from ©-© and (JT5J) by the identity 



<91nE 
dint 



dlnE 
dint 



cHnE 
dhiy 



(17) 



Along with hn(r) in (JSJ and (JTUJ) we require also 
the pair correlation functions hn (r)\= h,2i(r)] between 
solute-solvent pairs and h 2 2(r) between pairs of solutes, 
again in the limit of infinite dilution. The results are 
given here, with the detailed derivations in Appendix A. 

As shown in the appendix, 

^eJ) (r>l), (r>2), (18) 



hu(r) 



h\2{r) 



where D = hi 2 {l) / h\\{V) is a constant, independent of 
r. As also shown in the appendix, 



/ll 2 (l) 



1 



l/(c°a) 



(19) 



with c and a as defined in ([T} , ([3)) , and (0} ; while from © 
and ITUl). 



"■""=<;h£ 



i 



(20) 



Prom (1191) and (|20l) we then have the constant _D. In the 
meantime, h 22 (l) contains the solute-solute interaction 
energy, —£22- We define 



(21) 



analagously to the definition of a in (1161) . Then as further 
shown in the appendix, 



Ml) = -l + [l + /in(l)]. 



6-1 / 1 + ca 



1 + c a a 



(22) 



We also have /ii 2 (0) = h 2 2(°) = -1. From ([9]), ([ID]), 
and (|18 p - l|2"2"|) we now have all three correlation functions 
hn(r), hi 2 (r), h 22 {r) for all r. 

These results were derived both by the evaluation 
of the appropriate constrained partition functions and 
by transfer-matrix methods. The two derivations are 
equally lengthy; the appendix outlines the derivation via 
transfer matrices. 

Like ft-n(r) in ^ and (fTU[) . here, too, h\ 2 {r) and h,22{r) 
are independent of Z when the thermodynamic state of 
the system is specified by the temperature t and solvent 
density y, and so are the same as for the linear chain, 
Z = 2. 

The analogs of the Ornstein-Zernike relation for a two- 
component fluid are the Kirkwood-Buff relations 
In the limit in which the solute 2 is infinitely dilute in 
solvent 1, this requires for the present model that 



tx/n 



00 

v 2 /v + N r h 12 (r), 

r=Q 



(23) 



where 1*2 is the partial molecular volume of the solute 
at infinite dilution in the solvent. This is the analog of 
the Ornstein-Zernike relation (fl2)) with hn{r) in place 
of h\i(r) and ^2/^0 m place of V\/vq = 1/y. This i^/^o 
is in turn related by thermodynamic identity to the rate 
of ch ang e of solubility E with the density y of the sol- 
vent [18j |. 



— = -t- 

VQ 7T 



91nE 
dlny 



- 1 



so from 



it is required that 



00 

r=0 



91nE 
dhiy 



(24) 



(25) 



One may verify from (fT2J), ([15]), and ([18 ]l -([20 ]) . with 
^n(O) = hi 2 (0) = — 1, that the required ([2"5"]) is satis- 
fied. This is an important consistency test. 

From ([T8]) we have h 2 2{r) / h\ 2 (r) = 
hi 2 {r)/hn(r) (r > 2). This may be recognized to 
be a special case of a general principle definiti vely 
discussed in an important paper by Evans et al. [19|, 
who show that a similar relation holds asymptotically as 
r — >• 00 for any mixture, not only in the limit in which 
some of the components are present at infinite dilution. 
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In the present context it holds not only for r — >• oo 
but for all r > 2. One sees also in (fT5|) that hx2(f) 
and h22{f) have the same exponential decay length £ 
as in the solvent, given in (|10l) . That is because the 
correlations between a given solvent molecule and a 
molecule of an infinitely dilute solute, or between two 
such solute molecules, propagate through the solvent 
and so must have the range of the solvent-solvent 
correlations themselves. 

In the low-solubility states that are our primary in- 
terest, the common ratio D in (fig)) , evaluated as D = 
h\2{l)/hi\{l) via p^|) and (|2"0|) . proves to be very large 
(and negative), so the amplitude of /122M is then much 
greater than that of h\2(r), which in turn is much greater 
than that of hn(r). As a consequence, the amplitude 
of 7122(7") is enormously greater than that of hn(r). 
This property was earlier observed in a one-dimensional 
model of molecules interacting with square-well poten- 
tials [18j, [20( . It is likely to be generally true in solutions 
of solutes of low solubility, as a reflection of the effective, 



solvent-mediated attraction between them. 

The effective attraction between solute molecules is 
manifested also in the second osmotic virial coefficient, 
B. In a continuum model, with solute species 2 at infi- 
nite dilution, this would be given by —1/2 of the integral 
of /122 (V) over all space. In the present lattice model this 
is 

00 

- = --^^22(0 (26) 

with N r in ([TTj) . With Ji22(r) as above, and from the 
earlier formulas, the summation may be done explicitly 
and the result expressed as 

B = B +B 1 +B 2 , (27) 

where these are the contributions from r = 0, r = 1, and 
r > 2, respectively: 



Bo 
Bx 

B2 

vo 



1 

2' 



— Z(Z - l)D 2 
2 V ' 



i-1 

y 



1 

- - 1 

y 



Q-l 
Q + l 

Q-l 
Q + l 



„b-l 



1 



COL 



1 + c a a 
1 



i-(*-i)£r 



(28a) 
(28b) 

(28c) 



In the states of interest here, B/uo will prove to be large 
and negative, with the major contribution coming from 
Bi/vq; i.e., from the solute pair in contact. If we think 
of the negative B as arising from the "dimerization" [2l1 — 

m, 



so 



2 solute monomers 



solute dimer, 



with "monomer" and "dimer" both dilute in the solvent, 
then — B is the "equilibrium constant" for this dimeriza- 
tion. 

Let AG P , AS P , and AH p be the net changes in free 
energy, entropy, and enthalpy when a solute molecule is 
transferred at fixed p and T from an ideal gas phase where 
its number density is p2(i-g-) to the solution where its 
number density is p 2 (soln.). We define the constant- 
pressure solvation free energy and entropy AG* and 
AS; by AG* = AG P - fcTln [p 2 (soln.)/p 2 (i -gj] and 
AS* = AS p + kln [p 2 (soln.)/p 2 (i.g.)], while AH p itself 
is the constant-pressure solvation enthalpy ■ These are 
related by AGJ; = AH P -TAS;. When the ideal gas is in 
equilibrium with the solution, p 2 (i . g . ) is the same as the 
common activity Z2 in the two phases, while AG P = 0, 



ag; 



-feTlnE 



(29) 



with the solubility S as defined in (1131) . In the limit of 
infinite dilution, S is independent of the separate val- 
ues of P2 and Z2\ it depends on the nature of the solute 
and solvent and on their interactions but is otherwise a 
function only of the thermodynamic state of the solvent 
(Henry's law). The solvation entropy AS* and enthalpy 
AHp as defined above are then given by [7| 



as: 



AH P 
kT 



d_ 

df 



(TlnS) 



1. 



T 



ainS 
dT 



(30a) 
(30b) 



where 9 is the dimensionless coefficient of thermal expan- 
sion defined in ([5]) and the —1 is the analogous quantity 
for the ideal gas. 

By (|29|) . AG*/fcT is a measure of how low the solu- 
bility of the solute is. As remarked in the Introduction, 
a low solubility will be reflected in a strongly attractive 
potential of mean force between solutes. We concentrate 
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FIG. 1: Coordination number Z from fits of ()31[) with right- 
hand side to simulation data of Paschek 271] for Xe in various 
water models as solvent. 



t/t c = | ln9 w 0.44. In the states of interest, the choice 
a < 1/2 proves to be essential for the solubility to be low. 
The reason for this will be seen clearly in the asymptotic 
limit derived in Section [IVI 

From the formulas in Section [ill the liquid at the co- 
existence curve at this t then has the properties 

y = 0.999268, tt = 0.000127 
£ = 0.414 

X = 6.16 x 10~ 7 , 6 = 6.09 x 10" 3 . (33) 

This is a very high, nearly close-packed density and a very 
low pressure. The exponential decay length is very short: 
between r and r + 1 the correlation functions h\\(r) and 
h 12 {r) for r > 1, and /122(f) for r > 2, decay by the factor 
exp(— 1/£) = 0.089. The compressibility and coefficient 
of thermal expansion are both very low, as in real liquids. 

From the formulas in the present section, with values 
of the parameters as above, 



on the solvent-mediated part, W(r), of this potential; 
i.e., on the part remaining after the direct solute-solute 
interaction is subtracted. In particular, we wish to relate 
AG* to W(l), the contact value of W(r), for we expect 
these to be closely correlated when the solubility is low, 
as found in an earlier lattice model of hydrophobic ef- 
fects Q and confirmed in simulations of realistic models 
of aqueous solutions of non-polar solutes [27J. Indeed, 
as shown in Appendix B, we find these in the present 
Bethe-lattice model to be related by 



-w(i) + ag; = pv , 



(31) 



where, from (j2T)]) and (|22l) . but now without the factor 



-W(l)/kT _ 



1 



- 1 



1 



Q + i 



1 



ca 



1 + c a a 

'(32) 

In the states of most interest here, pvQ is very much 
smaller than either AG* or -(Z/2)W(1), so the right- 
hand side of (|3"Tj) may often be replaced by 0. When 
Paschek's simulation data [27| for Xe in SPC, SPC/E, 
TIP3P, TIP4P, and TIP5P water, in the temperature 
range 275 to 375K, are fit by (|3 1 1) with right-hand side 0, 
most of the resulting Z lie between 4.2 and 5.8 (Fig. fl}. 
Such values of the coordination number plausibly reflect 
the water structure. 

We now illustrate numerically the formulas in this and 
the preceding section by evaluating them in a represen- 
tative state in which the solute is of low solubility in the 
high-density liquid solvent that is in coexistence with its 
low-density vapor at a low temperature. We recall the 
abbreviations in ((T|), (fT6|) . and (f2~Tj). We choose for il- 
lustration Z = 3, a = 1/3, b = 1/9, and t = 1/5. 
(The relation b = a 2 is "Berthelot's rule"). With this 
Z the critical t from ([5]) is t c = 1/ In 9, so we have chosen 



E = 0.0865, (dmE/dlni)^ = 2.74 
D = -269.2 
B/v = -12.9 
{B Q /v = 0.5, Bi/vq = -11.9, B 2 /v = -1.6) 
AG*/fcT( = -lnE) = 2.447, W(l)/kT = -1.631 



ASl/k = -0.698, AH p /kT = 1.749. 



(34) 



This, as anticipated because of the high solvent density 
and low temperature (and because a < 1/2), is a low sol- 
ubility E, in the range of solubilities of small non-polar 
solutes in water. Also as anticipated, reflecting the low 
solubility, is a large negative D, which is the common ra- 
tio of h,22(r) to /112(f) for r > 2 and of h\2{r) to h\\{r) for 
r > 1. As a consequence, h22(r)/hn(r) = D 2 w 72,000. 
So great a ratio of amplitudes should be readily dis- 
cernible by experiment or simulation if one could mea- 
sure the correlations at large enough r. For example, with 
these values of the model's parameters, hn (2) = 6 x 10 -6 
while /i22(2) = 0.4; so, despite /122 and hu having the 
same formal decay length, /122 remains substantial long 
after hu has become negligible. That D is negative 
reflects the effective repulsion between solute and sol- 
vent associated with the low solubility, while the solvent- 
solvent interactions and the effective, solvent-mediated 
solute-solute interactions are both attractive. The sec- 
ond osmotic virial coefficient is large and negative, again 
reflecting the effective attraction between solutes that 
accompany the low solubility. As remarked earlier, the 
main contributor to the large negative B is the effective 
(direct and solvent-mediated) solute-solute attraction at 
r = 1. We note also that with the present Z = 3 the 
relation between W(l) and AG* required by the iden- 
tity (|3"Tj) is satisfied with pvo/kT (= ir/t) negligibly small 
compared with both -\ZW{l)/kT and AG*/kT. 

As a final remark, we observe that (dT./dt)^ > 0, 
AS** < 0, and AH p > 0. In these respects the low sol- 
ubility here then contrasts with that in the hydrophobic 
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effect, where (<9E/<9T) P < 0, AS; < 0, and AH P < 0. 
There, the sofubifity decreases with increasing T, and 
the sofubifity is iow because the sofvation entropy is suf- 
ficiently unfavorable to overweigh a favorable solvation 
enthalpy, whereas in the present model the solubility in- 
creases with increasing T, and the solubility is low be- 
cause the entropy and enthalpy of sofvation are both un- 
favorable. In the present model, too, there are states in 
which (9E/<9T) p < 0, but they are not the low-solubility 
states of interest. 



IV. MEAN-FIELD AND LOW-TEMPERATURE 
LIMITS 

The formulas in Sections [TT] and Mil simplify greatfy in 
the mean-field and low-temperature limits. 

The mean-field limit is that in which Z — > oo and 
en — > with the product Zen = <fi fixed. In this limit, 
from t in JT]) and from ([5]), 



= 4fcT c , 



and then from ©-©, 



pv 



fcTln- 



To l-y 



2y 2 



T c l-y 

1 In-r- 

_± i-v 



(35) 

(36) 
(37) 

(38) 



The equations of the coexistence curve and spinodal, 
from © and ([7]), are 



T 
T 



(coex . ) 



4y(l — y) (spinodal) 



(39) 
(40) 



and from (|15j) the solubility becomes 



(i - y)e 



4ayT a /T 



(41) 



Equation (|4ip may be understood from the potential- 
distribution formula. The factor (1 — y) is the prob- 
ability that the test solute molecule is in a previ- 
ously empty cell, Zy in the mean-field limit is the 
average number of neighboring cells occupied by sol- 
vent molecules, and the mean-field limit of the condi- 
tional average (exp(— ^ /kT)) is then exp(Zyei2/kT) = 
exp(Zyaen/kT) = exp(ya0/fcT) = exp(4yaT c /T). We 
see in this approximation that E always decreases with 
increasing T at fixed y (although not at fixed p). 

From (|2"9")l and (|3~Tj) . in this mean-field limit, 



Z 



1 



1 



W(l) = ~y(y-2a)h AG; = kTln-^-ay<j>, 
z 2 " l-y 

(42) 

which, with (|36p . is seen to satisfy the general (|3"TT) . 

A more relevant and more interesting limit is that of 
low temperature, which is equivalently that of high den- 
sity if we continue to specify the thermodynamic state 
of the liquid solvent as that at its coexistence curve, 
where it is in equilibrium with a low-density vapor. Then 
from ^ , asymptotically, as t and 1 — y both approach 
on the coexistence curve, they are related by 



(43) 



Now from the formulas of Section HlH in this limit, with 
Z > 2 (otherwise there is no coexistence curve) and a < 1 
(the 1,2 attraction weaker than the 1,1 attraction), 



and with b 



E 
D 

kT 



-Z{±-a)/t 
e (Z/2-a)/t 

Z W(1) Z{\-a) 
2 kT t 



(44) 
(45) 

(46) 



Bo 



Bi 



±Ze^l\ 



Eh 

r 



— Z(Z-l)e 
2 v ' 



[-iZ+2(l-o)]/i 



(47) 



From (j44|). it is now clear why a < 1/2 was a condition 
for low solubility in the low-temperature, high-density 
solvent states of interest here. Also, since Z > 2, we 
see from (|45[) why D is negative and of large magnitude, 
thus making the amplitude of h 2 2(r) hugely greater, by 
the factor D 2 , than that of hn(r) for all r > 2. We note 
that the large — D is correlated with the low solubility 
via -D ~ %-(z-2a)/Z{i-2a)^ w fa c ^ wit h a < 1/2, is a 

negative power of E. We see from (|47| why Bi/vq is neg- 



ative and of large magnitude. From (fiST) . we see that (|3"Tj) 
holds with the pvo on the right-hand side negligible com- 
pared with AG; and -\ZW(1) in this limit. [Cf. the 
dimensionless tt{— pvo/^ii) in (|33|)] . Numerically, with 
Z = 3, a = 1/3, t = 1/5, these asymptotic formulas yield 
the approximations y « 0.99945, E 0.082, D -341, 

AG* /kT |W(l)/fcT w 2.5, B x /v Q « -13.8, and 

B2/V0 « —1.3, which may be compared with their exact 
values in (EH) and (l34l). 
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V. SUMMARY 

We have studied a lattice model of a solution in which a 
solute is at infinite dilution in a liquid solvent. The lattice 
is a Bethe lattice (Cayley tree) of coordination number 
Z . The model parameters are ultimately chosen to make 
the solute solubility S, defined in (|13p . low, which then 
results in a strong, solvent-mediated attraction between 
solute molecules. 

The properties of the solvent alone, which are required 
in the later parts of the paper, are summarized in Sec- 
tion [TTJ The main new results are the properties of the 
model solution obtained in Section Mil and in the ap- 
pendices, and the low-temperature asymptotic limits in 
Section El 

The solute-solute, solvent-solute, and solvent-solvent 
pair-correlation functions hiiif), h\ 2 {r), and hn(r) 
are found in (TT51) to be related by h 2 2(r) / 'h\2{r) = 
hi2{r)/hn(r) for all r > 2. This is remarked to be a 
special case of an important statistical-mechanical the- 
orem on the asymptotic (r — > oo) behavior of the pair 
correlations in mixtures, as discussed by Evans et al. [19j. 
In the low-solubility, high-solvent-density states of inter- 
est, the common value D of those two ratios is found 
to be negative and large in absolute value. The result 
is that the amplitude of /122(f) is orders of magnitude 
greater than that of hn(r), and so is still discernible at 
distances beyond those at which hu{r) has become neg- 
ligible. This is in agreement with what was observed in 
an earlier model [ia . l2(ij |. This strong effective attrac- 
tion between pairs of solute molecules reflects the low 
solubility; i.e., the positive (thus, unfavorable) constant- 
pressure solvation free energy AG* = — /cTlnE. This is 
seen also in the second osmotic virial coefficient B, which 
in the states of interest is large and negative. The main 
contributor to B in these states is h 2 2(l), the solute- 
solute pair correlation function at contact, which shows 
the propensity to dimerization. 

The effective attraction between solutes that accom- 
panies low solubility is seen also in the quantitative re- 
lation pip , which relates the solvation free energy AG* 
to W(l), the solvent-mediated part of the potential of 
mean force between pairs of solute molecules at contact. 
The pressure term pvo on the right-hand side is usually 
negligible compared with both AG* and — \ZW(V) and 
may thus often be taken to be 0. Near linearity of the re- 
lation of AG* to W(l) was observed in an earlier lattice 
model as well Q , and confirmed in computer simulations 
of aqueous solutions of non-polar solutes [27j • Fitted val- 
ues of the coordination number Z are then in the range 
4.2 to 5.8. 

In the states of primary interest here, the solubility 
increases with increasing temperature at fixed pressure. 
The entropy of solvation AS* is negative and the en- 
thalpy AH p is positive, so AS* and AH p both contribute 
to the low solubility. These are in contrast to the hy- 
drophobic effect as seen in aqueous solutions of non-polar 
solutes, where £ decreases (AG* increases) with increas- 



ing temperature, and AH p < 0, so the solubility is low 
only because AS* is sufficiently negative to overweigh it. 
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Appendix A: Pair correlation functions 



Here we derive (0, dHJ) and ([221), with a solvent 1 at 
density y and a solute 2 at infinite dilution. We define 
P m > (r) as the probability that if a lattice site is occupied 
by a molecule of species 77, a site r lattice steps from it is 
occupied by a molecule of species r/'. "Species" means 
unocuppied. For a Cayley tree, with a solute molecule at 
a specified site, the probability P 2 i(r) of finding a solvent 
molecule at a site r steps away from it satisfies 

P 2 i(r) =P 2 i(r-l)Pii(l) + P 2 o(r-l)P i(l)- (Al) 

In addition, the P VV '(r) satisfy P21 (r— l)+P2o(r — 1) = 1, 
and the solvent density y at any site is related to the 
density at its neighbouring sites through: 

y = yP 11 (l) + (l-y)P i(l). (A2) 

P m ' (r) is related to h m i (r) by 

P 7iri ,{r)=y[l + h m ,{r)}, (A3) 

so from (|AT]) and (M\i , 

h2i(r)=d r - 1 h 21 (l), (A4) 

with h 2 i(l)[= h 12 {l)} from (19]) and d = P u (l) - P i(l). 
These also hold when solute and solvent are identical, 



MO 



%!(!), 



from (|A4l) and (jA5l) , 



M*0 = Ml) _ n 
h n (r) Ml) ~ ' 



(A5) 



(A6) 



which is the first of the identities f(T8j) . 

Next we calculate the pair distribution function 
522 {r) [= 1 + h<22 (r)J, which is given by the potential dis- 
tribution theorem Q: 



922 (r) 



-P[UL B (.0) + U'(r)] 



(A7) 



where /3 = 1/kT, 6ij is the Kronecker delta, and B (k) 
is the total interaction energy between a solute molecule 
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FIG. 2: A branch of the Cayley tree with coordination number 
Z — 3 and a solute on its base site. 



and the surrounding fluid. The denominator appearing 
in (gll is related to jUJ through £ = (e^^sM). T h e 
numerator of dSI, (e^ IC, « M+E, ^ (r)1 ), will be calcu- 
lated using transfer matrices [Hj]. To simplify the repre- 
sentation of the matrices' elements we will follow ref. [29| 
and use the so called normalized partial partition func- 
tions x n , which correspond to the ratio between partition 
functions of branches with a molecule 77 and a hole on its 
base site (Fig. These functions are related to the 
quantities previously defined in the text by 



X2 = 



II 



{l-y)a 
1 + ac a 



1 



(A8a) 
(A8b) 



The transfer matrices connecting different lattice types 
are given by: 



V 



(11) 



= [V] 



r)r)' 

C W x (Z-2)( n +r,' ■)/2 z (n+ V ')/2 



~v (21) ~ 




'v {12) 




1)1)' 





(A9) 



= Sor,^ 1 -^ 'x[ z - 2W/2 z*/ 2 , (A10) 

where r\ and 7/ are occupation variables restricted to un- 
ocuppied sites (0) and solvent (1), and where z is the 
solvent activity, which satisfies zxf =y/(l — y). 

To obtain an expression for g 22 (?") in terms of a trace, 
we connect sites and r with ghost links whose elements 
are: 



A (H) 
A (22) 



>/>/' 



= x ftl +, l'V 2 z (V+i)')/2 

■ - A A JZ-i)(2-i)-i)') 



With the above definitions, and for r > 2, 
1 lT:{A^V^V r - 2 V {1 ^} 



922 (r) = 



E 2 



lr{A( n )V r } 



(Alia) 
(Allb) 

(A12) 



After some algebra, and by using (|A8a[) and (lA8b[) . one 
finds: 



322 (r) = 1 + 



1 



y 



D 2 e~ r ^., 



(A13) 



where £ is given by (jlOl) and 

,2 _ Mi) 



Mi) ! 



(A14) 



which completes (|T8|) . 

Since eq. (|A13|) is valid only for r > 2, it is necessary 
to calculate /i22(l)- We recalculate (|A7|) for r — 1 using 



(A15) 



-^, B (Q) + ^(i)]\ = JrjA^V^} 
Tr{A( n )^} : 



with V (22) defined as 

"y(22) 



'/'/' 



(A16) 

After some manipulations (| A7|) . at contact, becomes 

Together with eqs. © and (fT5j) . this completes the eval- 
uation of /ii2(?") and h 22 {r). 



Appendix B: The relation between W(l) and AGp 

The relevant part of the potential of mean force be- 
tween two solutes does not contain the solute-solute in- 
teraction, being exclusively the solvent-mediated part of 
the effective interaction. At contact, r = 1, this is given 
by 



W(l) = -kTln[g 22 (l)] + e 22 . 
From equation (|A17I) it becomes: 

'{l+y)x 1 + y 



W{1) = 2kTlnx 2 + kTln 



Xi 



(Bl) 



(B2) 



The standard solvation free energy of the solute 2, from 
flU]), P9|- and (lA8bll . is: 



AG* = -A:Tln(l - y) - Z/cTlna; 2 , 



(l+y)xi +y 

X\ 



(B3) 



so from eqs. (|B2I) and (IB3 



-W(l) + AG; = fcTln- 



1-1/ 



(B4) 

Using eqs. (HJ)- (J3J) and (|A8al) one may show that the 
right hand side of (|B4[) is pvg , which then gives (|3ip . 



9 



G. S. Rushbrooke, Introduction to Statistical Mechanics 
(Oxford University Press, 1949), pp. 300-304. 
T. L. Hill, Statistical Mechanics (McGraw-Hill, 1956), 
pp. 348-353. 

G. M. Bell and D. A. Lavis, J. Phys. A: Gen. Phys. 3, 
427 (1970). 

G. M. Bell and D. A. Lavis, J. Phys. A: Gen. Phys. 3, 
568 (1970). 

N. A. M. Besseling and J. Lyklema, J. Phys. Chem. B 
101, 7604 (1997). 

C. D. Eads, J. Phys. Chem. B 106, 12282 (2002). 

B. Widom, P. Bhimalapuram, and K. Koga, Phys. Chem. 
Chem. Phys. 5, 3085 (2003). 

N. Guisoni and V. B. Henriques, J. Phys. Chem. B 110, 
17188 (2006). 

T. L. Beck, M. E. Paulaitis, and L. R. Pratt, The Poten- 
tial Distribution Theorem and Models of Molecular Solu- 
tions (Cambridge University Press, 2006), Chap. 7. 

D. M. Rogers and T. L. Beck, J. Chem. Phys. 129, 
134505 (2008). 

M. A. A. Barbosa and V. B. Henriques, Phys. Rev. E 77, 
051204 (2008). 

C. Buzano and M. Pretti, J. Chem. Phys. 119, 3791 
(2003). 

J. C. Wheeler and B. Widom, J. Chem. Phys. 52, 5334 
(1970). 

J. G. Kirkwood and F. P. Buff, J. Chem. Phys. 19, 774 
(1951). 



[15] B. H. Zimm, J. Chem. Phys. 21, 934 (1953). 

[16] A. Ben-Nairn, Water and Aqueous Solutions: Introduc- 
tion to a Molecular Theory (Plenum Press, 1974), p. 142. 

[17] A. Ben-Nairn, J. Chem. Phys. 67, 4884 (1977). 

[18] P. De Gregorio and B. Widom, J. Phys. Chem. C 111, 
16060 (2007). 

[19] R. Evans, R. J. F. Leote de Carvalho, J. R. Henderson, 
and D. C. Hoyle, J. Chem. Phys. 100, 591 (1994). 

[20] P. De Gregorio, J. C. Toledo, and B. Widom, Mol. Phys. 
106, 419 (2008). 

[21] E. E. Tucker and S. D. Christian, J. Phys. Chem. 83, 426 

(1979) . 

[22] P. J. Rossky and H. L. Friedman, J. Phys. Chem. 84, 587 

(1980) . 

[23] C. Chipot, P. A. Kollman, and D. A. Pearlman, J. Com- 
putational Chem. 17, 1112 (1996). 

[24] H. Liu and E. Ruckenstein, J. Phys. Chem. B 102, 1005 
(1998). 

[25] O. Coskuner and U. K. Deiters, Z. Phys. Chem. 220, 349 

(2006) . 

[26] O. Coskuner and U. K. Deiters, Z. Phys. Chem. 221, 785 

(2007) . 

[27] D. Paschek, J. Chem. Phys. 120, 6674 (2004). 

[28] S. R. A. Salinas, Introduction to statistical physics 

(Springer- Verlag, New York, 2001), Chap. 13. 
[29] C.-K. Hu and N. S. Izmailian, Phys. Rev. E 58, 1644 

(1998). 



